Conserving diversity in Irish plant–pollinator networks

Abstract Beneficial insects provide valuable services upon which we rely, including pollination. Pollinator conservation is a global priority, and a significant concern in Ireland, where over half of extant bee species have declined significantly in recent decades. As flower‐visiting insects rely on flowering plants, one way to conserve and promote pollinator populations is to protect high‐quality habitat. We analyzed the structure of insect–flower interactions from multiple habitat categories in a large database of interactions from Ireland. Our primary goals were to compare spatial and temporal variation in Irish network structures, compare Irish networks to published networks from other countries, and provide evidence‐based recommendations for pollinator conservation in Ireland by identifying well‐visited plant species that may promote high pollinator diversity, abundance, and functional complementarity. Habitat types within Ireland differed substantially: seminatural grasslands had the highest pollinator species richness and largest number of unique pollinator species, while intensively managed habitats exhibited negative asymmetry (more plant than pollinator species). This negative asymmetry is notable because most plant–pollinator networks exhibit a positive asymmetry. Within intensively managed habitats, agricultural and urban habitats differed. Urban habitats had the highest number of non‐native plant species while agricultural habitats had the lowest pollinator species richness. We also found Irish networks varied across the growing season, where July had the highest plant and insect species richness. When comparing Irish networks to published networks from other countries, we found Irish networks had a higher ratio of plant species to pollinator species, and that this difference was most evident in agricultural habitats. This ratio means the typical network asymmetry (more pollinator than plant species) was flipped (more plant than pollinator species) in the Irish network. We conclude that conserving seminatural grasslands in Ireland will be an essential component of pollinator conservation and identify thirty‐five plant species important for restoring seminatural habitats.

variation in Irish network structures, compare Irish networks to published networks from other countries, and provide evidence-based recommendations for pollinator conservation in Ireland by identifying well-visited plant species that may promote high pollinator diversity, abundance, and functional complementarity. Habitat types within Ireland differed substantially: seminatural grasslands had the highest pollinator species richness and largest number of unique pollinator species, while intensively managed habitats exhibited negative asymmetry (more plant than pollinator species).
This negative asymmetry is notable because most plant-pollinator networks exhibit a positive asymmetry. Within intensively managed habitats, agricultural and urban habitats differed. Urban habitats had the highest number of non-native plant species while agricultural habitats had the lowest pollinator species richness. We also found Irish networks varied across the growing season, where July had the highest plant and insect species richness. When comparing Irish networks to published networks from other countries, we found Irish networks had a higher ratio of plant species to pollinator species, and that this difference was most evident in agricultural habitats.
This ratio means the typical network asymmetry (more pollinator than plant species) was flipped (more plant than pollinator species) in the Irish network. We conclude that conserving seminatural grasslands in Ireland will be an essential component of pollinator conservation and identify thirty-five plant species important for restoring seminatural habitats.

K E Y W O R D S
agricultural intensification, conservation, flower-visiting insects, grasslands, plant-pollinator networks

| INTRODUC TI ON
Globally, ecological systems have been substantially altered by human activities (Ellis & Ramankutty, 2008). For example, the primary effect of agricultural food production at a local scale is the loss of species through the application of agrochemicals and conversion of land toward high-yield monocultures (Flohre et al., 2011;Grab et al., 2019;Matson et al., 1997;Tscharntke et al., 2005). Other forms of human activities alter ecosystems differently, for example, through the introduction and promotion of non-native ornamental species, such as in urban or developed landscapes (Godefroid & Koedam, 2007). Regardless of mechanism, habitat degradation has unfavorable implications for the preservation of beneficial biodiversity, particularly those species that provide ecosystem services which maintain sustainable human production systems (Biesmeijer et al., 2006;Hallmann et al., 2017;Lister & Garcia, 2018;Ollerton et al., 2014;Powney et al., 2019).
Complex mutualistic networks of plant-pollinator interactions provide valuable ecosystem services for agricultural production (Dainese et al., 2019;Klein et al., 2007;Rader et al., 2016) alongside their natural role in plant reproduction more broadly (Ollerton et al., 2011). The importance of pollinator diversity for agricultural yield has been demonstrated by many studies (Calderone, 2012;Dainese et al., 2019;Garibaldi et al., 2016). At the same time, agricultural land-use can have serious negative impacts on pollinator diversity . The impacts of human activities can be observed in highly coevolved networks of interacting plants and pollinators, for example, by showing increases in non-native plant and insect taxa over time (Mathiasson & Rehan, 2020). As species in such networks are mutually dependent, species loss can have cascading implications for the persistence of the species with which they interact. Changing environmental quality can also be reflected in structural changes, for example, by altering the asymmetry of insect-flower networks (Soares et al., 2017). Asymmetry is considered a fundamental attribute of mutualistic network structure, and may reflect important coevolutionary patterns (Bascompte et al., 2006), so changes to asymmetry may impact the stability of these interactions over time. The introduction of non-native species may also increase species richness and add novel interactions in a way that leaves the overall network structure unchanged, while changes in interactions between native species are not detectable at the network level (Lopezaraiza-Mikel et al., 2007).
It is thus essential to understand and mitigate the impacts of human land-use on beneficial pollinating insects. It is clear that the highest priority conservation action is protecting remaining intact, high-quality habitat. A secondary objective is improving the quality of degraded systems, for example, by identifying central organisms in plant-pollinator networks that support both a diversity and abundance of partners. As organisms contribute to community structure on multiple levels, it can be important to evaluate them using various metrics.
To better understand the state of plant-pollinator interactions within Ireland, we addressed three main aims. To understand how habitat quality affects plant-pollinator interactions in Ireland, we first compared networks within Ireland to one another. We compared the structure of networks in intensively managed (agricultural and urban), coastal dune, grassland (seminatural and managed), and woodland/shrubland habitats in Ireland. We further compared changes in these Irish networks across the growing season. Then, to view Irish networks in a broader context, we compared Irish network structures to published networks from other countries. Finally, we evaluated the role of individual species within the Irish networks to provide evidence-based recommendations for habitat restoration or conservation-oriented planning.  (Larkin & Stanley, 2021;O'Rourke et al., 2014;Power & Stout, 2011;Stanley & Stout, 2013  Agroecology, Biodiversity ecology, Community ecology, Conservation ecology, Landscape ecology insects for later identification in the laboratory when field identification was not possible.

| Irish data
These studies comprised more than 940 h of surveys of plantpollinator interactions across 709,721 m 2 and 119 study sites in four broad habitat categories (Table S1): (1) cultivated and built land (agricultural crops [including Miscanthus, winter wheat, oil seed rape]), silage and dairy pastures, and urban and suburban gardens and flower beds; (2) coastland (primarily coastal dune systems); (3) grassland and marsh (primarily seminatural and managed grasslands); and (4) woodland and scrub (primarily hedgerows and non-native shrubs). Our habitat classification was based on the Fossitt (2000) guide to Irish habitats, which has three separate levels of categorization. At the first level, there are 11 habitat types defined within Ireland, from which we sampled four thoroughly. It is worth noting that at the first level of categorization, managed (improved for agriculture) grasslands and seminatural grasslands are grouped together in the broad category of "grassland." We therefore further separated the habitats into a level 2 classification. There are 30 subcategories of habitat types at level 2, of which we sampled from 18 ( Figure S1; Table S1). For the level 2 classification, we were able to separate managed grasslands from seminatural grasslands. We selected the most extensive habitat types that represented the majority of the land surface in Ireland. A rarefaction analysis showed our sample coverage for these four habitats varies from 96.6% (coastland) to 99.3% (cultivated and built land) ( Figure S2). When comparing habitat types, we excluded records where the surveyor included more than one habitat type in a single survey.
From these datasets, we automatically excluded any data where the host plant or insect visitor were not resolved to at least genus.
Although most specimens were resolved to the species level, we included some specimens that had only been identified to genus, as some plant and insect taxa were difficult/impossible to resolve to species, often listed as "aggregates" (e.g., Taraxacum agg. and Bombus lucorum agg.). These studies excluded plants that were present but not visited by the target insects.

| Statistical software packages
We conducted all of our statistical analyses using R (R Core Team, 2013). To conduct a rarefaction analysis, we used the package iNEXT (Chao et al., 2014). To analyze the attributes of the bipartite networks, we used the package bipartite (Dormann et al., 2009), specifically calculating the weighted and unweighted degrees of plant species and insect visitor species, network level connectance, nestedness (using NODF), and asymmetry. In addition to the attributes calculated using the bipartite package, we conducted a functional complementarity analysis (Devoto et al., 2012) with the package vegan (Oksanen, 2011) and calculated centrality by creating a one mode, unweighted projection of plant and insect species in the networks using the package igraph (Csardi & Nepusz, 2006). To calculate the effect of floral display size on visitation, we used the package lme4 (Bates et al., 2014). We also used this package to determine whether visitation to non-native plants was different than to native plants, given their floral display. NMDS plots of flower-visitor communities in different habitat types and months of the year were constructed using the package vegan (Oksanen, 2011).

| Habitat type and phenology in Irish networks
We compared the network structure among four primary habitat types where interactions were recorded in Ireland, intensively managed (agricultural or developed), coastal dunes, grasslands (including both seminatural and managed), and woodlands/shrublands.
We quantified the species richness of plants and pollinators in each habitat, along with the number of non-native plant species. We constructed networks weighted by the abundance of insect visitors on plant species for each habitat type based on the interactions recorded and calculated plant and insect visitor species richness, and weighted and unweighted degree, as well as network connectance, nestedness (NODF), and asymmetry (Table S2). We then identified insect species common to all habitat types (Table S3), as well as those unique to each habitat type in this step (i.e., recorded in only one habitat type) (Table S4). We compared the flower-visitor communities in these habitats using a nonmetric multidimensional scaling (NMDS) ordination. To determine the relative role of agriculture versus other processes, we further divided habitats into (a) primarily agriculture versus urban intensively managed habitat types, and (b) seminatural versus managed (species-poor and used for agricultural purposes) grasslands, and evaluated their attributes. These subcategories were assigned according to level two of the Fossitt habitat classifications (Table S1).
Next, using the same data, we constructed separate networks for each of the five months which comprised most of the surveys (May-September). For this comparison, we constructed networks from data pooled across years for each month and reported the same measures calculated above (Table S2). For both habitat and phenology (i.e., month) data, we constructed NMDS plots to visualize any overlap in flower-visitor community structure.

| Null models for specialization metrics
To test the degree of network-level and species-level specialization, we compared the existing networks to null network models constructed using the nullmodel function in the package bipartite.
We used this function to generate 1000 null models for each habitat (the full Irish network, human-modified habitats, coastal dunes, grasslands, and woodlands) and month (the full Irish network, May, June, July, August, and September). Null networks are generated using the vaznull method, which randomizes the total number of interactions in the original network, while keeping the connectance constant (Vazquez & Aizen, 2003). We then calculated the mean network-level specialization (H 2 '), and specialization of the two separate groups of organisms (insects and plants) (d') across all null models (Blüthgen et al., 2006). We used z-scores to compare the true values to the null model averages.

| Comparing the Irish network to other plantvisitor networks
To determine whether the attributes of the Irish network were similar to other published flower-visitor networks, we downloaded publicly available network data from the Interactions Web Database and DataDryad (Table 3). We specifically downloaded networks of quantitative plant-insect visitor interactions with more than a total of 50 species (including both plant and insect species) represented, to ensure the networks were similar in size to our network, as some network attributes are sensitive to networks smaller than 50 species (Dormann et al., 2009). We found nine such networks in the Interaction Web Database and one in DataDryad.
As some of these networks included non-syrphid dipterans, we repeated the analysis while excluding these groups from the datasets, but saw no significant difference in our results ( Figure S6; Table S5).
Using the downloaded datasets, we calculated the species richness of plants and insect visitors, and the following network measures: weighted and unweighted average degrees for plant and insect species, as well as the connectance, nestedness (NODF), and asymmetry of the networks. We selected these measures because they are commonly evaluated in flower-visitor network studies and have been well-studied, with theoretical underpinnings relating to stability, co-evolution, and robustness to species loss (Table S2). Moreover, they have been suggested as conservation targets and potentially an important component of monitoring (Tylianakis et al., 2010). As we were comparing a single value (i.e., the value of the Irish network), we determined that the Irish network was significantly different when its value was more than two standard deviations from the mean of the other published networks.

| Ranking Irish plant and insect species
There tends to be an inherent background relationship between flower-visitor abundance and the size of the floral display (Russo et al., 2013). Flower-visitors may exhibit a preference for plant species by visiting them in greater abundance than expected based on their availability, described here as the size of the floral display. In contrast, a plant species with a large floral display but low visitor abundance may be considered less preferred. Therefore, using generalized linear mixed effects models (GLMMs), we tested the relationship between floral area and the abundance and richness of insect visitors to determine whether more insects were attracted to a larger floral display. We used a log-log transformation to evaluate the relationship between abundance, richness, and floral area. To account for repeated measures, site and date of sampling were used as random effects in these models. We also investigated whether visitation to a plant family was greater than anticipated given the background abundance and floral display of the species within that family.
The contributions of individual plant and insect species to the structure of the Irish networks were ranked according to several different metrics (Table S2). We categorized these species by functional complementarity (calculated with a hierarchical clustering algorithm) (Devoto et al., 2012), average abundance of visitors (weighted degree) and average partner species richness (unweighted degree), visitation rate (abundance of visitors divided by the floral display), and node longevity (duration of species activity in the network across the weeks of the summer), as well as centrality in unweighted one-mode projections of the network (closeness and betweenness centrality) (Russo et al., 2013). We normalized these values for all species on a scale of 0-1, where 0 represents minimum values and 1 represents maximum values (i.e., feature scaling), to give different measures equal weight. For plant species, we also created subsets of the network data and recalculated the network parameters separately based on visitation by the following insect subgroups: bees, syrphid flies, and Lepidoptera. Depending on management objectives, it may be preferable to optimize different measures to focus on plant species conservation and restoration. Land managers cannot construct insect communities in the same fashion as plant communities, but this information may be useful for targeting conservation objectives toward more central pollinators, or for protecting key agricultural pollinators.
We used Pearson Correlation Coefficients to test for correlations between the measures listed above. Most of these measures significantly correlated with one another, with the exception of visitation rate, which could be therefore considered an independent measure (Table S6). Functional complementarity was negatively correlated with all of the other measures, and likely also provides a contrasting measure of a species' role in the community.
To test whether non-native plant species differed significantly from native plant species in terms of their role in the networks, we generated GLMMs containing the following fixed effects: the status of the plant (native or non-native) and the responses of functional complementarity, weighted and unweighted degree, node longevity, visitation rate, and betweenness and closeness centrality. To correct for multiple testing, we used a Bonferroni correction (alpha = 0.007).
For these models, site and week of the year were inserted as random effects to account for repeated measures, and the fixed effects were categorical (native or non-native species and habitat type). We used a Poisson distribution in these models.

| Habitat type and phenology in Irish networks
Although intensively managed habitats had the highest plant species richness, they also had the highest proportion of non-native plant species. Moreover, they had proportionally lower insect species richness, leading to the most negative asymmetry; intensively managed habitats had 2.44 times as many flowering plant species as insect visitor species (Table 1). The only habitat where the asymmetry was not negative (i.e., more visitor species than plant species) was coastal dune systems, which also had the highest connectance. In wood and shrubland habitats, there remained nearly twice as many plant as insect species (Table 1). The grassland habitats (including both seminatural and managed) had the highest insect species richness overall, along with the greatest average abundance and species richness of insect visitors per plant species. Due to their high plant species richness, and corresponding low insect species richness, intensively managed habitats also had the highest abundance of insects per visitor species, and plant species richness per visitor, along with the highest nestedness (Table 1).
There was substantial overlap in insect visitor species across habitats; the most abundant visitors were found in all habitats and visited the largest number of plant species (Figure 1, Table S3). On the other hand, only four plant species were found in all four habitat types (Rubus fruticosus agg. L., Senecio jacobaea L., Taraxacum agg. L., and Trifolium repens L.), each of which has a global distribution, suggesting wide ecological tolerance. The highest proportion of unique insect visitor and plant species was found in grasslands (Table S4), which also had the highest visitor species richness (Figure 1), average visitor abundance (effect size 0.58, z value 2.36, p = .02), specialization (H 2 '), and specialization asymmetry (d'), relative to the null model comparison (Figure 2a,b). Approximately 40% of both insect and plant species recorded in grasslands were only found in that habitat type (Table S4). The grasslands contrasted with the intensively managed habitats, which had the lowest proportion of unique species (Table S4). The woodland/shrubland habitats had the lowest specialization and the full network had the most negative specialization asymmetry (Figure 2a,b). We also found that woodland/  We also include subsets of the main habitat types where they are distinguished by human land-use. For example, the intensively managed habitats can be divided into agriculture only or urban habitats, while the grasslands can be separated into seminatural or amenity (improved for nonagricultural usage).
The highest value for each measure is indicated in bold.
As the proportion of non-native plant species increased, asymmetry and nestedness tended to decrease across the habitat types (Table 1; Figure S3). However, agricultural habitat types were distinct; they had a low proportion of non-native plant species, a strongly negative asymmetry, and a high nestedness (Table 1; Figure S3). Grasslands managed for agricultural use through the application of fertilizers and herbicides had a negative asymmetry and high nestedness, and the highest proportion of non-native species of the habitat subsets (Table 1; Figure S3).
We compared network measures for different months of the flowering season and found that the month of July stood out as having the highest abundance-related measures and total richness, but that the month of August had a higher average partner species richness and connectance. The highest nestedness and the most positive asymmetry were found in the month of September, although that month also had the lowest species richness overall.

| Comparing Irish networks to other networks
Relative to the other networks analyzed here, including those of comparable latitude, species richness, size, and habitat types, and

F I G U R E 2
The z-scores of the full Irish network (black) and habitat subsets as compared to null models for the mean specialization index (H 2 ') of the network of the habitats (a, b) and months (d, e) within the Irish network as compared to null model simulations. The grasslands have the highest overall specialization across the full network and least negative specialization asymmetry (relative to the mean of the null models). NMDS ordination plots of the flower-visitor communities in the different habitats (c) and months (f).  The highest value for each measure is indicated in bold. (Table S2): insect species richness (a), plant species richness (b), insect abundance per insect species (c), insect abundance per plant species (d), plant species richness per visitor (e), visitor species richness per plant (f), connectance (g), nestedness (h), and asymmetry (i). Each plot contains the median (dark line), the upper and lower quartiles (box), the range (dotted lines), and various outliers (dots). The full Irish network (black) and four habitat types (intensively managed [primarily agriculture] = orange, coastal dunes = yellow, grasslands = green, and wood and shrublands = blue) are indicated by circles when they do not differ from the other published networks and stars when they differ significantly (by more than two standard deviations from the mean).      (Table 1). Thus, the difference between the grasslands overall and other networks were attributable to managed grasslands.

| Relative species importance
Overall, the strongest predictor of insect visitation in the Irish plantpollinator network was the size of the floral display (floral abundance × inflorescence size) (Figure 4) We ranked plant species in the network by eight measures, and then recorded the top ten species according to each metric. Thirtyfive plant species were ranked in the top ten of at least one measure ( Figure 5). Four of these (11%) were non-native, which was substantially lower than their prevalence in the overall network ( (Figure 5f). The non-native Fuchsia cultivar had the highest visitation rate, followed by Sambucus nigra L., Cirsium vulgare (Savi) Ten., C. arvense, and L. saxatilis (Figure 5a).
Among the distinct insect groups, R. fruticosus agg. supported the greatest abundance of syrphid and bee visitors, but Lotus corniculatus L. supported the greatest abundance of lepidopterans. Bees, moths, and butterflies also preferred the plant family Fabaceae over its background prevalence ( Figure 6).
When ranking the insects, twenty-eight insect species were ranked in the top ten of all metrics of "importance," 53.6% of which were syrphid flies, 42.9% bees, and 3.6% Lepidoptera ( Figure S5).
While the proportion of syrphids in the top ten reflected their frequency in the database (54.7% syrphids), there were proportionally more bee species, but fewer Lepidoptera, ranked as highly important compared with their background prevalence (30.7% bee and 14.6% Lepidoptera).

F I G U R E 4
Scatterplot showing the significant relationship between log transformed floral area and log transformed average abundance of visitors in the full Irish network, with non-native plant species (blue) significantly less attractive than native species (red). Gray points on the graph (NA) represent plant species where the origin is unknown, either because they are cultivars or because they are not resolved to the species level. The gray shaded area around the lines indicates the 95% confidence interval of a linear relationship.
When ranking the insect visitors, Bombus pascuorum Scopoli, ranked the highest for every measure except for functional complementarity, whereas Halictus Latreille ranked the highest ( Figure S5).

B. pascuorum visited 127 of the 239 plant species (53%).
Non-native plant species represented 50% of plant species richness on average across surveys, but received only 11% of total visits.
They also had a significantly lower average flower-visitor abundance (effect size −0.83, t value −6.91, p < <.001, Figure 4). We found nonnative plant species differed significantly from native plant species in terms of their network roles (Table 4). They had a lower abundance and richness of visitors, bloomed for a shorter duration of the flowering season, and had a lower betweenness (but not closeness) centrality. Non-native species did not differ from native species in their functional complementarity, visitation rate, or closeness centrality (Table 4).

| DISCUSS ION
As the intensity and extent of human land-use increases (Tilman et al., 2011), it is important to consider how we can best protect communities of beneficial insects. More than half of Irish bee species have experienced significant declines in recent decades (Fitzpatrick et al., 2006). Clearly, the first step is to protect natural and seminatural habitats which retain beneficial biodiversity. In Irish landscapes, seminatural grassland habitats retained the highest proportion of unique species of both insect visitors and plant species: 40% of the visitor species in the grasslands were only found in that habitat and included four bee species listed as "vulnerable" by the recent red list of Irish bees (Andrena angustior Kirby, A. semilaevis Pérez, Bombus ruderarius Müller, and Lasioglossum nitidiusculum Curtis) (Fitzpatrick et al., 2006). There were only four plant species common to all habitat F I G U R E 5 Dotplots showing the thirty-five plant species in the Irish network that rank in the top ten of the different measures (visitation rate (a), average abundance of visitors (weighted degree, (b)), average species richness of visitors (unweighted degree, (c)), betweenness centrality (d), duration of bloom (node longevity, (e)), and functional complementarity (f)) ranked from highest (top) to lowest (bottom). Our analysis also illustrated patterns in the phenology of the network, showing the highest abundance and species richness of both plant and visitor species in July, while the network was more connected in August. September had the highest nestedness and asymmetry of the months, but also was represented by relatively few plant and insect species. A study in Pennsylvania (USA), showed a similar pattern of the highest species richness in the middle of the summer, but also found the highest nestedness at that time, while connectance was highest early in the season (Russo et al., 2013).
Interestingly, the full (year-round) network was more negatively asymmetric than any month alone, because the plants tended to be active for shorter durations of the summer than their insect partners. This overlap in visitor communities among the months can also be illustrated by the overlap in their NMDS ordination.
Most flower-visitor networks are highly asymmetric, with a much larger number of insect visitors than plant hosts. This asymmetry has implications for coevolutionary processes (Bascompte et al., 2006;Ramírez et al., 2011;Russo et al., 2019). In contrast, the compiled F I G U R E 6 A heat map of the Irish network depicting the relative preferences of different insect visitor groups (all together (Full), bee, hoverfly (Syrphid), and moth and butterfly (Lep)). These demonstrate the prevalence of members of each plant family in the top ten rankings relative to the background floral abundance of plant species in each plant family in the overall network. Blue represents a group that is underrepresented in the top ten rankings relative to its background abundance, and red highlights groups that are overrepresented relative to their background abundance. The values are relative proportions different from background abundance. For example, syrphids prefer species of the family Asteraceae at 20% higher than the background abundance of these plant species in the network. Note: For these models, we set the network measure as the response variable, and the status of the plant (non-native or native) as the fixed effect. The family for the models was Gaussian. We report the effect size, t value (Wald Statistic) and p value. To account for repeated tests, we use a Bonferroni correct alpha = 0.007. Significant effects at this alpha are highlighted and in bold.

TA B L E 4 Results of generalized linear mixed effects models comparing network roles of native and non-native plant species in the Irish network
data for the Irish flower-visitor network showed a strong asymmetry in the opposite direction: across all habitats, there were far more plant than visitor species. These differences were largely driven by intensively managed, especially agricultural, habitats in Ireland.
Managed grasslands had the highest proportion of non-native plant species, but did not have proportionally greater numbers of pollinator species. In fact, the managed grasslands are often deliberately seeded with agricultural species to increase their value for grazing or silage (Fossitt, 2000). Primarily agricultural habitats had few nonnative plants and relatively low pollinator species richness, suggesting non-native plant species were not contributing to the network structure in that habitat. Woodland/shrubland habitats in Ireland also had negative asymmetry, possibly because the habitats in this category were dominated by samples from hedgerows, which tend to line agricultural fields, and non-native shrubs, including ornamental plants found in nonagricultural intensively managed habitats.
Seminatural grasslands had high flower-visitor species richness and low proportions of non-native plant species, leading to a positive asymmetry more in line with what was observed in networks outside Ireland.
The Irish network was also more nested than other networks, possibly because of the highly asymmetrical interactions, as nestedness has been linked to asymmetry (Bascompte et al., 2003).
Moreover, it has been argued that nestedness relates to community persistence in mutualistic networks. On the other hand, the Irish network was significantly more nested than other networks we observed, and significant variation in nestedness among mutualistic networks calls into question its importance as a predictor of biodiversity. Moreover, nestedness in seed dispersal networks was higher in species-rich networks that were located in areas of high human impact (Sebastián-González et al., 2015) and may result from high environmental variability (Song et al., 2017).
Where habitat has been degraded, it is important to identify central plant species that support a large diversity and abundance of insects and those that support a unique complement of insect  (Pywell et al., 2005).
The relative species richness of Syrphidae in the Irish network suggested they may be undervalued contributors to pollination services in Ireland. However, in terms of sheer abundance and number of plant species visited, the more abundant and generalized bees (especially B. pascuorum and Bombus lucorum agg.) outranked Syrphidae and thus were more "important" by our measures.
Similarly, to the plant species rankings, this is driven by the broad distribution and high abundance of these bumble bee species. These two common bumble bees were responsible for 30% of the total observed visits. In Ireland, highly generalized bumble bee species represent both the core of the network structure and also the vast majority of interactions.
Flower-visiting insects provide a valuable ecosystem service and though pollination was not directly measured by this study, other studies have demonstrated the importance of flower-visiting insect diversity to supply pollination services in agricultural systems (Dainese et al., 2019;Garibaldi et al., 2013). If there is a detrimental net effect of agriculture on flower-visitor diversity, there are serious implications for the sustainability of pollination services for pollinator-dependent crops and terrestrial ecosystems. Ireland has a long history of agriculture (Brown, 2007), and the highest proportion of land in Europe dedicated to agricultural production (Eurostat Agricultural Census, 2010), but the impact of agriculture is not restricted to Ireland. As the human population and demand on agricultural production increase, both the intensity and extent of agriculture can only increase globally (Tilman et al., 2011). Ireland may exemplify the future of plant-pollinator interactions globally.

ACK N OWLED G M ENTS
Identifications for insect specimens were aided by the National

CO N FLI C T O F I NTE R E S T
The authors have no conflict of interest.

DATA AVA I L A B I L I T Y S TAT E M E N T
The data for this manuscript are available at DataDryad at the following DOI: https://doi.org/10.5061/dryad.kwh70 rz47.